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Abstract 

This  paper  presents  a  detailed  numerical  investigation  of  the  transport  and  electrochemical  phenomena  involved  during  the  operation  of  a 
single  proton  exchange  membrane  (PEM)  fuel  cell  on  reformate  feed  with  a  view  to  developing  optimal  design  and  operating  conditions.  A 
one-dimensional  non-isothermal  model,  validated  with  experimental  data,  is  utilized  to  evaluate  the  fuel  cell  performance  over  a  wide  range 
of  design  and  operating  parameters  that  affect  the  thermal  response  and  water  management.  Based  on  a  systematic  parametric  analysis  on  the 
various  physical  and  electrochemical  phenomena,  feasible  operating  regimes  and  optimal  design  conditions  are  identified  with  the  objective 
of  maximizing  the  power  density  subject  to  constraints.  Overall,  this  paper  illustrates  a  methodology  for  using  physics-based  models  for  cell 
design  and  optimization. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Proton  exchange  membrane  (PEM)  fuel  cells  are  promis¬ 
ing  power-generation  sources  for  mobile  and  stationary  ap¬ 
plications.  In  a  typical  PEM  fuel  cell,  a  polymer  membrane  is 
inserted  between  an  anode  and  a  cathode  electrode  to  form  a 
membrane  electrode  assembly  (MEA),  which  is  further  sand¬ 
wiched  between  two  bipolar  plates  housing  the  flow  channel, 
as  shown  in  Fig.  1 .  A  thin  catalyst  layer  exists  between  each 
of  electrodes  and  the  membrane,  referred  to  as  the  anode  and 
cathode  catalyst  layer,  respectively.  The  advantages  of  PEM 
fuel  cells  include  the  ability  to  provide  high  current  densities 
at  relatively  low  operating  temperature  and  pressure,  quick 
start-up,  and  pollution-free  operation  [1].  The  high  cost  and 
low  reliability  of  fuel  cells,  however,  are  the  limiting  factors 
for  their  widespread  commercialization.  A  good  understand¬ 
ing  of  the  effects  of  the  design  and  operating  conditions  on 
the  voltage  losses  is  required  to  reduce  the  capital  cost  and 
improve  the  reliability.  Accordingly,  two  major  groups  of 
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parameters  may  be  identified,  namely  the  operating  parame¬ 
ters,  e.g.,  temperature,  pressure,  reactant  stoichiometry,  and 
gas  composition,  and  design  parameters,  e.g.,  thickness  of 
cell  components,  electrode  porosity,  and  platinum  catalyst 
loading. 

Water  and  thermal  management  are  critical  to  the  overall 
cell  performance  [2-7].  To  maintain  ionic  conductivity,  the 
membrane  in  a  PEM  fuel  cell  requires  adequate  humidifi¬ 
cation,  which  raises  the  critical  issue  of  water  management. 
During  practical  operation  of  fuel  cells,  both  the  gas  streams 
are  humidified  to  ensure  the  proper  membrane  hydration. 
However,  excessive  water  will  accumulate  in  the  electrode 
pores  and  result  in  the  electrode  flooding,  which  degrades  the 
cell  performance  by  preventing  the  reactants  from  reaching 
the  catalyst  sites.  In  general,  higher  operating  temperature  is 
desirable  due  to  decreased  mass  transport  limitations  and  in¬ 
creased  electrochemical  reaction  rates;  at  the  same  time,  high 
temperatures  may  lead  to  increased  mass  transport  losses  due 
to  the  increase  in  water  vapor  pressure.  A  careful  design  of 
the  cell  and  its  operating  parameters  is  therefore  imperative 
to  balance  such  competing  constraints. 

In  the  last  decade,  efforts  on  modeling  and  computer  sim¬ 
ulation  of  PEM  fuel  cells  have  been  directed  towards  bet- 
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Nomenclature 

a  active  reaction  area  per  unit  volume 

(cm2  cm-3) 

bfcO  back-to-forward  CO  adsorption  ratio  (atm) 
bfh  back-to-forward  FL  adsorption  ratio  (atm) 

Cj  molar  concentration  of  species  i  (mol  cm-3) 

Cpj  specific  heat  at  constant  pressure  for  species  i 
(J  mol-1  K-1) 

Dj  diffusivity  of  species  i  in  membrane  (cm2  s-1) 

Dij  diffusivity  of  species  i  in  species  j  (cm2  s-1) 

ee  electrode  porosity 

/fCe] i  cell  potential  difference  (V) 

F  Faradays  constant,  96487  (Cmol-1) 

/ivap  enthalpy  of  vaporization  (J  g- 1 ) 
z'm  catalyst  layer  membrane  phase  current  density 
(A  cm-2) 

4  catalyst  layer  solid-phase  current  density 

(A  cm-2) 

z'o  exchange  current  density  (A  cm-2) 

I  operating  current  density  (A  cm-2) 

j(X)  total  reaction  rate  (A  cm-3) 

ji(X)  reaction  rate  of  species  i  (A  cm-3) 
kfc  electrode  forward  CO  adsorption  rate  times  2 F 
(A  cm-2  atm-1) 

kfl,  electrode  forward  FL  adsorption  rate  times  2 F 
(A  cm-2  atm- 1 ) 

kp  hydraulic  permeability  (cm2) 

L  characteristic  length  scale  (cm) 

m  catalyst  loading  (gem-2) 

Nj  molar  flux  of  species  i  (mol  cm-2  s-1) 
p  pressure  (atm) 

pa  anode  pressure  (atm) 

R  universal  gas  constant  8.314  (J  mol-1  K-1) 
AS  entropy  change  (J  g-1  K-1) 

Sli  Sherwood  number 

/ci  catalyst  thickness  (cm) 

/ei  electrode  thickness  (cm) 

/me  membrane  thickness  (cm) 

T  temperature  (K) 

Wi  molar  mass  of  species  i  (gmol-1) 
x  position  coordinate  along  the  cell  thickness 

Xi  mole  fraction  of  species  i 

Greek  symbols 

aa  apparent  transfer  coefficient  for  anodic  reac¬ 
tion 

ac  apparent  transfer  coefficient  for  cathodic  reac¬ 

tion 

8  membrane  expansion  coefficient 

i]d  electro-osmotic  drag  coefficient 

k  thermal  conductivity  (W cm~ 1  K~ 1 ) 

X  membrane  hydration  coefficient  (moles  of  wa¬ 

ter  per  mole  charge  sites) 


/x  viscosity  (g  cm- 1  s- 1 ) 

cf)m  membrane  phase  potential  (V) 
solid-phase  potential  (V) 
a  electronic  conductivity  (£2-1  cm-1) 

6\  hydrogen  coverage 

04  CO  coverage 

tOj  mass  source  for  species  i  (gem-3) 

Subscripts 

an  anode 

ca  cathode 

cat  catalyst 

el  electrode 

8  gas 

l  liquid 

me  membrane 

ref  reference 

Superscripts 

eff  effective 

i  boundary  i 


ter  understanding  of  the  cell  operation  and  performance. 
Bernardi  and  Verbrugge  [8,9]  and  Springer  et  al.  [10,11] 
studied  the  steady-state  isothermal  operation  of  PEM  fuel 
cells  using  a  one-dimensional  model,  and  assuming  perfect 
membrane  hydration.  Baschuk  and  Li  [12]  studied  the  ef¬ 
fects  of  variant  degree  of  water  flooding  in  the  cathode  elec¬ 
trode/catalyst  layer  on  the  overall  cell  performance  using  a 
steady-state  one -dimensional  approach.  Their  study  demon¬ 
strated  that  flooding  has  a  significant  impact  on  cell  perfor¬ 
mance.  Rowe  and  Li  [13]  performed  a  complete  study  on  the 
water  and  thermal  management  on  PEM  fuel  cells  using  a 
steady-state  one-dimensional  approach,  and  prescribed  wa¬ 
ter  vapor  mole  flux  at  the  interfaces  of  electrodes  and  catalyst 
layers. 

A  two-dimensional  model  of  transport  phenomena  in  PEM 
fuel  cells  was  presented  by  Gurau  et  al.  [14],  and  the  effects 
of  oxygen  and  water  vapor  mole  fraction  distribution  on  the 
cell  performance  were  illustrated.  Um  et  al.  [15]  proposed 
a  transient,  single -phase  two-dimensional  model  for  electro¬ 
chemical  and  transport  processes  in  a  PEM  fuel  cell.  The  heat 
management  was  neglected  by  assuming  isothermal  opera¬ 
tion.  Amphlett  et  al.  [16,17]  studied  the  transient  response 
of  a  fuel  cell  stack  by  performing  a  global  heat  and  mass 
balance  analysis,  while  ignoring  the  details  of  electrochemi¬ 
cal  phenomena  inside  the  cell.  A  two-dimensional  isothermal 
model  was  developed  by  Wang  et  al.  [18]  to  study  the  two- 
phase  flow  and  transport  in  the  cathode  side.  The  liquid  wa¬ 
ter  formation  and  its  effects  on  the  reaction  rate  and  species 
transport  were  analyzed.  You  and  Liu  [19]  and  You  [20] 
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Fig.  1.  Schematic  illustration  of  a  PEM  fuel  cell. 


presented  a  two-dimensional  steady-state  isothermal  model 
which  couples  the  flow  species,  potential  and  current  density 
distributions  in  an  individual  fuel  cell.  Although  the  two- 
phase  flow  was  restricted  to  the  cathode  electrode,  the  study 
shows  that  the  incorporation  of  two-phase  flow  in  the  model  is 
essential. 

While  modeling  of  the  individual  transport  and  electro¬ 
chemical  phenomena  has  been  reported  collectively  in  the 
literature  and  provides  insight  on  fuel  cell  operation,  a  sys¬ 
tematic  study  to  quantitatively  determine  the  optimal  condi¬ 
tions  from  physics-based  models  is  lacking.  In  view  of  this, 
a  methodology  for  model-based  design  and  optimization  is 
presented  for  the  first  time  based  on  the  systematic  paramet¬ 
ric  studies  and  consideration  of  a  few  illustrative  constraints 
on  cell  operation.  The  model  adopted  in  the  present  study 
combines  the  transport  model  from  Rowe  and  Li  [13]  and  the 
CO  poisoning  kinetics  from  Springer  et  al.  [21],  as  briefly 
presented  in  the  next  section.  The  numerical  model,  vali¬ 
dated  with  the  available  experimental  and  numerical  results 
in  the  literature,  is  used  to  conduct  a  parametric  exploration 
in  terms  of  the  operating  and  design  parameters.  Based  on 
the  parametric  studies  and  specified  constraints  on  the  maxi¬ 
mum  allowable  temperature  difference,  maximum  allowable 
cell  voltage  and  minimum  desirable  membrane  hydration, 
operating  windows  are  identified  on  the  current  density  as  a 
function  of  the  various  parameters. 


2.  Mathematical  model 

To  introduce  the  methodology  for  model-based  design  and 
optimization  of  a  PEM  fuel  cell,  a  physical  model  for  the 
transport  and  electrochemical  phenomena  is  needed.  In  this 


section,  the  transport  model  from  Rowe  and  Li  [13]  and  the 
electrochemical  kinetics  model  from  Springer  et  al.  [21]  are 
combined  to  predict  the  performance  of  a  PEM  fuel  cell  oper¬ 
ating  on  reformate  feed.  The  model  considers  the  five  distinct 
regions  between  the  flow  channels,  as  shown  schematically 
in  Fig.  1 .  The  cell  is  considered  to  be  operating  at  steady  state, 
and  since  the  primary  aim  of  the  study  is  to  present  a  method¬ 
ology  for  design  and  optimization  of  cells  using  physics- 
based  models,  the  discussion  is  limited  to  a  one -dimensional 
modeling  in  the  direction  along  the  cell  thickness.  The  mod¬ 
eling  further  includes  the  effects  of  carbon  monoxide  (CO) 
poisoning  of  the  catalysts,  as  is  prevalent  in  fuel  cells  oper¬ 
ating  on  a  reformate  feed,  but  is  often  neglected  in  fuel  cell 
modeling.  The  governing  equations  for  the  electrode,  cata¬ 
lyst  and  membrane  regions  are  presented  in  the  following 
sections. 

2.1.  Electrodes 

Fuel  cell  electrodes  are  typically  made  of  porous  carbon 
paper  or  cloth,  which  serves  to  transfer  the  reactant  species 
and  to  conduct  electrical  current.  The  mathematical  model 
is  obtained  by  considering  the  conservation  of  the  species, 
momentum  and  energy.  Since  the  viscous  force  and  pres¬ 
sure  gradient  are  assumed  to  be  zero  in  the  electrode  regions, 
the  solution  of  the  one-dimensional  steady-state  momentum 
equation  yields  constant  velocity  of  the  feeding  gas  mixtures. 
Following  the  development  by  Rowe  and  Li  [13],  the  govern¬ 
ing  equations  may  be  written  as 

d  Ni  a>i 

Species  :  - =  —  (1) 

d.r  Wj 
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Energy  : 


— A-eff 


d  2T 

dx2 


NiWiCpJ  +  NiW3cP'i 


d  T 

x  — 
dx 


1 2 

+  /lvap^3  ^  (^) 


where  the  subscript  i  denotes  the  ideal  gas  species  i,  Nj  the 
molar  flux  in  the  x-direction  in  Fig.  1 ,  o>,  the  mass  source  term, 
Wj  the  molecular  weight,  kctt  the  effective  thermal  conductiv¬ 
ity,  T  the  temperature,  cpj  the  specific  heat  at  constant  pres¬ 
sure,  /V|  the  molar  flux  of  liquid  water,  h  vap  the  enthalpy  of 
vaporization  for  water,  /  the  current  density,  <rett  the  effective 
electrical  conductivity  and  the  summation  is  performed 
with  respect  to  all  the  gas  species  in  the  mixture. 

In  this  study,  the  feed  streams  are  considered  to  be  refor¬ 
mate  fuel  at  the  anode  side  and  humidified  air  at  the  cathode 
side.  Hence,  the  gas  species  i  are  defined  as  1  =  O2,  2  =  N2, 
3  =  HoOtg)  (i.e.,  water  vapor)  for  the  cathode  electrode;  and 
1  =  H2,  2  =  CO2,  3  =  H20(g)  and  4  =  CO  for  the  anode 
electrode.  Note  that  the  inclusion  of  the  CO  species  in  the 
modeling  forms  the  primary  distinction  of  the  present  model 
relative  to  that  reported  by  Rowe  and  Li  [13].  Since  no  elec¬ 
trochemical  reaction  occurs  in  the  electrode  regions,  the  mass 
source  term  &>,-  is  nonzero  only  for  the  water  vapor  species 
i  =  3.  Evidently,  Eq.  (1)  indicates  that  the  molar  flux  for  the 
species  O2,  N2,  H2,  CO2  and  CO  are  constants.  To  solve  the 
water  vapor  flux  N3,  the  mass  source  term  &>3  (correspond¬ 
ing  to  the  vaporization/condensation  of  the  water  species) 
is  determined  in  terms  of  the  temperature  and  species  con¬ 
centration  X3,  which,  in  turn,  is  determined  by  the  Stefan- 
Maxwell  equation  [13].  Furthermore,  the  electrical  potential 
in  the  electrode  in  the  electrode  solid,  tPs,  is  calculated  by  the 
Ohm’s  law,  and  the  unknowns  in  the  electrode  regions  are  T, 
N3,  X3  and  <PS. 


2.2.  Catalyst  layers 

Catalyst  layers  are  considered  to  be  a  mixture  of  mem¬ 
brane,  platinum  catalyst  (solid)  and  void  space  in  this  study. 
Electrochemical  reactions  in  the  catalyst  regions  are  coupled 
with  the  transport  of  mass  and  energy,  resulting  in  a  potential 
gradient  across  the  cell.  The  overall  half-cell  electrochemical 
reaction  in  the  cathode  catalyst  layer  may  be  written  as  [13]: 

02  +  4H+  +  4e“  — *  2H20(1)  (3) 

where  the  oxygen  species  and  the  electrons  from  the  cathode 
electrode  region  react  with  the  protons  from  the  membrane 
region  to  form  liquid  water.  For  a  PEM  fuel  cell  operating 
on  reformate  feed  at  the  anode,  the  electrochemical  reactions 
in  the  anode  catalyst  layer  involve  the  competing  adsorption 
processes  of  CO  and  H2,  and  are  represented  by  the  four 
processes  expressed  in  the  following  equations  [21]: 


(4) 

(5) 


2(M-H)  — *  2H+  +  2e“  +  2M  (6) 

H20  +  (M-CO)  — ►  M  +  C02  +  2H+  +  2e~  (7) 

The  bidirectional  arrow  in  Eq.  (4)  indicates  that  a  free  carbon 
monoxide  molecule  may  be  adsorbed  to  a  vacant  catalyst  site, 
M,  to  form  a  CO  molecule  in  the  adsorbed  state,  (M-CO);  in  a 
reverse  process,  the  adsorbed  CO  molecule  may  be  desorbed 
to  yield  a  free  CO  and  a  vacant  site.  Similarly,  a  hydrogen 
molecule  may  experience  a  reversible  chemisorption  process 
with  two  vacant  catalyst  sites  to  form  two  adsorbed  hydrogen 
atoms.  The  two  processes  in  Eqs.  (6)  and  (7)  generate  the 
current  corresponding  to  the  electrochemical  oxidation  of  the 
adsorbed  hydrogen  atoms  and  carbon  monoxide  molecule, 
respectively.  Clearly,  the  CO  and  H2  molecules  compete  with 
each  other  for  the  vacant  catalyst  sites,  and  high  concentration 
of  CO  in  the  fuel  may  prevent  the  adsorption  of  the  hydrogen 
molecules,  leading  to  the  so-called  CO  poisoning  effect. 

In  this  section,  the  subscript  notation  for  the  species  in¬ 
volved  in  the  anode  catalyst  layer  is  1  =  H2,  2  =  H+,  3  = 
H2Od),  4  =  CO,  5  =  CO2,  and  that  for  the  cathode  catalyst 
layer  is  1  =  O2,  2  =  H+,  3  =  H20(1).  The  governing  equa¬ 
tions  for  the  various  electrochemical  and  transport  processes 
are  derived  by  the  application  of  conservation  laws  for  the 
species  and  energy,  along  with  the  Bulter-Volmer  equation 
for  the  electrochemical  reactions,  the  Nernst-Planck  equa¬ 
tion  for  the  flux  of  aqueous  species  in  the  membrane,  and  the 
Ohm’s  law  for  electron  transfer  in  the  solid. 


2.2.1.  Anode  catalyst  layer 

The  conservation  equations  of  species  and  energy  in  the 
anode  catalyst  layer  account  for  the  electrochemical  reactions 


Species  H2  : 
Species  H+  : 
Species  H->0(1) 
Species  CO  : 
Species  C02  : 


Energy 


follows: 

diVi 

j]  (x) 

(8) 

dx 

IF 

diV2 

7t  W  +  74(x) 

(9) 

F — -  = 
dx 

dx 

d  N3 

Mx) 

(10) 

dx 

2  F 

diV4 

74(x) 

(11) 

dx 

2  F 

d/V5 

74  (x) 

(12) 

dx 

IF 

q 

dx2  4 

■(e 

NiCpyWi 

5  / 

dr 

71  +  74 

(T  AS) 

dx 

2  F 

a 

—  Ol  +  /4)(^>S  —  ^m)  +  7^ 


CO  +  M  (M-CO) 
H2  +  2M  2(M-H) 


(13) 
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where  the  reaction  rates  j\(x)  and  74  (x)  for  Ho  and  CO,  re¬ 
spectively,  may  be  given  by  the  Butler- Volmer  equation  as 


71 M  =  «*oef0i-^f 

ci 


exp 


/«a  F 

v  RT 


—  exp 


-acF 

RT 


f  C4 

74  to  =  ai^OA-^ 


exp 


aaF 
—  (^s  -  <Pm) 


—  exp 


-acF . 


RT 


tos  -  <t> m) 


(14) 


(15) 


The  parameter  ;m  is  the  catalyst  layer  membrane  phase  cur¬ 
rent  density,  which  is  related  to  the  proton  molar  flux,  N2,  via 
the  Faraday  constant,  F,  as  im  —  FN2 ■  In  the  energy  equation, 
Eq.  (13),  AS  is  the  entropy  change  for  the  cathode  reaction, 
(ps  and  0m  are  the  electrical  potential  in  the  catalyst  solid 
phase  and  catalyst  membrane  phase,  respectively,  and  /celf 
is  the  effective  electrical  conductivity  of  the  membrane.  The 
reaction  rates,  71  to  and  74 (x),  depend  on  the  catalyst  reac¬ 
tive  surface  area  per  unit  volume,  a,  the  reference  exchange 
current  density,  if^1,  at  the  reference  oxygen  concentration, 
cf,  the  fuel  concentrations,  c  1  and  C4,  and  the  transfer  co¬ 
efficients  aa  and  ac.  Depending  on  the  characteristics  of  the 
half-cell  reactions  and  the  material  properties  of  the  catalyst 
layers,  the  transfer  coefficients,  aa  and  ac,  take  on  distinct 
values,  each  in  a  range  between  0  and  1.  Following  Rowe 
and  Li  [13],  both  transfer  coefficients  are  taken  to  be  unity 
in  this  study,  which  yields  an  approximate  Tafel  slope  of  70 
mV  decade-1 . 

Since  both  CO  and  kb  are  oxidized  to  produce  H+,  the 
species  equation  for  proton  flux,  Eq.  (9),  includes  the  sum 
of  71  (x)  and  74(x).  The  species  consumption/generation  in 
the  electrochemical  reaction  of  carbon  monoxide  is  repre¬ 
sented  by  Eqs.  (10)— (12).  In  the  energy  equation,  Eq.  (13), 
the  summation  ^  includes  four  non-ionic  species,  and  the 
sum  7i(x)  +  72 (x)  is  also  introduced  for  the  overall  anode 
electrochemical  reaction  involving  Hi  and  CO.  Note  that  the 
reaction  rate  71  (x)  is  proportional  to  the  coverage  of  hydrogen 
molecules,  0 1,  which  is  defined  as  the  fraction  of  the  catalyst 
reactive  surface  area  covered  by  the  adsorbed  hydrogen.  Sim¬ 
ilarly,  the  reaction  rate  of  the  CO  species  is  considered  to  be 
proportional  to  the  CO  coverage,  64,  and  CO  concentration 
C4  in  Eq.  (15).  Considering  the  adsorption,  desorption  and  re¬ 
action  processes  in  Eqs.  (4)-(7),  the  coverage  9\  and  64  may 
be  obtained  from  a  kinetic  analysis  for  mass  balance  [21]: 


Table  1 

Kinetic  parameters  for  the  electrochemical  reactions  of  CO 


bfco  (atm) 

1.5  x  nr8 

bfh  (atm) 

0.5 

Fec  (A  cm-2) 

lxitr8 

Feh  (A  cm-2) 

4.0 

kfc  (A  cm-2  atm-1) 

10 

£fh  (A  cm-2  atm-1) 

4000 

8(AGCO)/RT 

7.8 

8(AEh)/RT 

4.6 

dissociative  adsorption  near  CO-occupied  catalyst  site.  Since 
the  species  Fb  and  CO  exhibit  different  adsorption/desorption 
kinetics,  with  the  CO  following  a  more  complex  Temkin  ki¬ 
netics,  the  values  of  S(AGco)  and  <5(A£’h)  are  different,  as 
seen  in  Table  1  [21], 

The  six  conservation  equations  for  the  anode  catalyst 
layer,  Eqs.  (8)— (13),  introduce  four  additional  unknowns, 
namely  ci,  C4,  0S  and  <£m,  which,  in  turn,  are  determined 
by  the  Nernst-Planck  equation  and  Ohm’s  law  [13].  Note 
that  10  unknowns  are  involved  in  the  anode  catalyst  region: 
N i,  ;m,  N$,  N4,  N5,  T,  ci,  C4,  0m  and  @s.  The  governing 
equations  for  the  cathode  catalyst  layer  follow  those  given 
by  Rowe  and  Li  [13],  and  are  not  repeated  here  for  brevity. 

2.3.  Membrane 


The  membrane  of  a  PEM  fuel  cell  is  generally  made  of  a 
perflourosulfonate  polymer,  which  acts  as  a  proton  conductor 
when  saturated  with  water.  Based  on  the  model  assumptions, 
only  two  species  are  present  in  the  membrane  region,  i.e.,  the 
liquid  water  and  the  proton.  Since  no  chemical  reaction  takes 
place  in  the  membrane,  the  species  conservation  indicates 
constant  fluxes  for  both  species.  The  flux  of  water  is  deter¬ 
mined  by  the  net  effect  of  electro-osmotic  drag,  diffusion  due 
to  concentration  gradient  and  convection  due  to  pressure  gra¬ 
dient.  The  flux  of  protons  is  described  by  the  Nernst-Planck 
equation,  which  is  further  rearranged  to  be  in  the  form  of  the 
membrane  potential  in  this  study.  The  conservation  equations 
of  energy  and  species  in  the  membrane  region  read  [13]: 

Energy  :  +  -jWqW 1)  =  ^  (18) 

dx-  dx  k 


Potential  : 


d0m  _  fm  s—(  3  \dx 

dx  k  F  V  1  +  <V, )  dx 

+  f(i)N,  (19) 


kfcX4 p(  1  -  04  -  01 )  -  0fc&fc04  -  74  =  0  (16) 

kfhX\p(l  -  04  -  0i )2  -  bfhkfhOj  —  71=0  (17) 

where  the  CO  adsorption-to-desorption  ratio,  b(c.  is  a  func¬ 
tion  of  the  free  energy  variation,  S(  AGcoX  f°r  the  CO  adsorp¬ 
tion/desorption  process,  and  Fb  adsorption  rate,  depends 
on  the  change  in  activation  energy,  S(AEn),  for  hydrogen 


H20  flux : 


dcj 

Ni  =  -Di-f-e 
dx 


mem 

w 


Vd  1 
F 

(20) 


where  the  subscript  1  denotes  the  liquid  water  species,  h  the 
specific  enthalpy,  k  the  proton  conductivity  of  the  membrane, 
S  the  membrane  expansion  coefficient,  X  the  membrane  hy- 
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dration,  D\  the  diffusion  coefficient  for  liquid  water  in  the 
membrane,  e™em  the  volume  fraction  of  water  in  the  mem¬ 
brane,  kp  the  hydraulic  permeability  of  the  membrane,  //  the 
viscosity  of  liquid  water  and  ijc\  the  electro-osmotic  drag  co¬ 
efficient.  Eqs.  ( 1 8)— (20)  are  solved  for  the  three  unknowns: 
T,  and  N\,  and  the  readers  are  referred  to  Rowe  and  Li 
[13]  for  the  values  of  various  physical  properties  involved. 

The  governing  equations  formulated  in  Sections  2. 1-2.3 
are  solved  with  the  following  boundary  conditions  specified 
at  the  inlets  of  the  anode  and  cathode  electrodes  (i.e.,  at  the 
points  a  and  f  in  Fig.  1):  temperature,  composition  of  the 
reactant  gas  mixtures,  pressure  and  flow  rate  in  terms  of 
stoichiometry.  It  must  be  pointed  out  that  the  water  vapor 
flux,  N%  could  be  calculated  by  considering  the  condensa¬ 
tion/vaporization  processes  in  the  porous  electrode  regions. 
However,  due  to  the  difficulty  in  solving  two-phase  flow  in 
the  porous  media,  the  values  of  N$  at  the  electrode/catalyst 
interfaces,  i.e.,  points  b  and  e,  are  set  to  be  10%  of  the  corre¬ 
sponding  flux  of  the  reactant  mixture,  as  suggested  by  Rowe 
and  Li  [13]  for  convenience.  An  algorithm  developed  by  Fan 
and  White  [22]  was  implemented  to  predict  the  polarization 
curve,  temperature  distribution,  species  concentration  and 
flux  under  various  operating  conditions.  The  physical  proper¬ 
ties  and  kinetic  data  adopted  in  the  simulations  are  available 
in  Refs.  [13,21], 


3.  Results  and  discussion 

The  validation  of  the  current  model  is  first  examined 
through  comparison  with  numerical  and  experimental  results 
in  the  literature.  Using  the  parameter  combination  for  the 
base  case  listed  in  Table  2,  the  predicted  polarization  curve 
is  compared  to  the  numerical  results  from  Rowe  and  Li  [13] 
in  Fig.  2  for  the  case  of  pure  hydrogen  at  3  atm  as  fuel.  In  the 
figure,  the  line  corresponds  to  the  present  simulation  results, 
while  the  discrete  circle  markers  denote  numerical  predic¬ 
tion  by  Rowe  and  Li.  The  results  of  the  model  in  this  paper 
are  seen  to  be  in  good  agreement  with  those  of  Rowe  and  Li 
for  the  range  of  current  density  considered.  The  cell  poten¬ 
tial  decreases  monotonically  with  increasing  current  density, 
owing  to  increased  ohmic  loss  in  the  membrane,  and  activa¬ 
tion  and  concentration  losses  in  the  catalyst  regions.  When 
the  current  density  increases  to  around  1 .0  A  cm-2,  the  con¬ 
centration  polarization  becomes  dominant  due  to  the  limita¬ 
tion  of  mass  transport  rate.  The  partial  pressure  of  oxygen 
at  the  catalyst  layer/cathode  interface  approaches  zero  when 
7=1.1  A  cm-2,  which  is  taken  as  the  limiting  current  den¬ 
sity  for  this  case. 

Fig.  2  also  compares  the  polarization  curve  prediction  with 
experimental  data  from  Springer  et  al.  [21]  for  a  PEM  fuel 
cell  operated  on  reformate  feed  with  three  CO  concentra¬ 
tions.  In  addition  to  Table  2,  the  kinetic  parameters  used  in 
the  calculation  for  the  electrochemical  reactions  of  CO  are 
given  in  Table  1 .  Again,  good  agreement  is  observed  between 
the  model  simulations  and  the  experimental  measurements. 


Table  2 

Base  case  model  parameters 

Cell  temperature,  Tc  (K)  353.0 

Anode  pressure,  pa  (atm)  3.0 

Cathode  pressure,  pc  (atm)  3.0 

Anode  stoichiometry,  £a  1.5 

Cathode  stoichiometry,  £c  3.0 

Relative  humidity,  RH  100 

Electrode  thickness,  te\  (|xm)  200 

Catalyst  layer  thickness,  tc\  (|Jim)  7 

Membrane  thickness,  tme  (|xm)  180 

Volume  fraction  of  membrane  in  catalyst  layer,  e^  0.45 

Volume  fraction  of  solid  in  catalyst  layer,  ef  0.5 

Anode  dry  gas  mole  fraction  (CO2/H2)  0.0 

Cathode  dry  gas  mole  fraction  (N2/O2)  3.76 

Reversible  potential,  Eoc  (V)  1.199 

Anode:  ai Qef  (A cm-3)  30000 

Cathode:  aiTff  (A cm-3)  0.0095 

Thermal  conductivity  electrode,  ke\  (W  cm-1  K-1)  0.026 

Thermal  conductivity  of  dry  membrane,  kme  (W  cm  1  K-1)  0.00673 

Thermal  conductivity  of  water,  £w  (W  cm-1  K-1)  0.0024 

(Sh  /e)/L  (cm-2)  20000 


Recall  that  the  CO  species  is  adsorbed  on  the  catalyst  sur¬ 
face,  blocking  the  electrochemical  reaction  of  the  hydrogen 
fuel.  Consequently,  smaller  cell  potential  is  seen  for  a  fixed 
value  of  current  density  when  CO  concentration  increases 
from  0  to  250  ppm  in  Fig.  2.  The  limiting  current  density  for 
the  three  cases  with  CO  poisoning  effects  reduces  to  around 
0.6  A  cm~2,  due  to  the  high  coverage  of  CO  species  and  low 
partial  pressure  of  the  hydrogen  at  the  catalyst  layer/anode 
interface.  Note  that  a  sharp  drop  of  7sceii  is  observed  between 
I  =  0.2  and  0.4A  cm-2  for  the  case  of  100  ppm  CO  concen¬ 
tration.  For  the  highest  CO  concentration  considered  in  this 
study,  i.e.,  250  ppm,  significant  loss  of  £cen  is  observed  as 
I  >  0.1.  Overall,  the  model  predictions  show  good  agreement 
with  the  numerical  and  experimental  results  in  the  literature. 
This  forms  a  reliable  basis  for  the  use  of  the  numerical  model 


Current  Density,  I  [A/cm2] 


Fig.  2.  Validation  of  the  present  model  with  numerical  result  of  Rowe  and  Li 
[13]  and  experimental  data  of  Springer  et  al.  [21]  for  the  CO  concentrations 
identified. 
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in  conducting  the  parametric  studies  discussed  in  the  remain¬ 
der  of  this  section. 

For  better  understanding  of  the  transport  and  electrochem¬ 
ical  processes  involved,  it  is  illustrative  to  first  examine  the 
distributions  of  the  species  concentration  and  flux,  tempera¬ 
ture  and  potential  across  the  fuel  cell  thickness.  To  this  end, 
Figs.  3  and  4  present  the  variables  of  interest  as  functions  of 
the  location  x  within  the  cell  at  three  current  densities.  The 
values  of  the  parameters  adopted  in  Fig.  3  are  listed  in  Table  2. 
Fig.  3(a)  shows  the  influence  of  current  density  on  the  con¬ 
centration  distribution  of  the  species  FF  and  O2.  For  the  case 
considered,  the  species  in  the  anode  electrode  are  the  hydro¬ 
gen  and  water  vapor,  which  corresponds  to  a  relatively  large 
value  of  diffusivity  for  FF  [13].  Thus,  the  Stefan-Maxwell 
equation  indicates  that  a  small  concentration  gradient  of  kb 
is  needed  in  the  porous  anode  to  maintain  the  hydrogen  flux 


to  the  anode  catalyst  layer.  However,  the  mutual  diffusivities 
of  the  oxygen  in  the  humidified  air  mixture  are  small,  lead¬ 
ing  to  relatively  large  c o2  gradient  in  the  cathode,  as  seen  in 
Fig.  3(a).  Note  that  the  positive  gradient  of  co  corresponds 
to  the  oxygen  flux  from  the  cathode  electrode  to  the  cath¬ 
ode  catalyst  layer.  Increase  in  the  current  density  calls  for  the 
increase  in  the  feed  stream  fluxes,  which,  in  turn,  is  respon¬ 
sible  for  the  increase  in  the  gradients  of  c o2  and  ch2  in  the 
electrodes  by  virtue  of  the  Stefan-Maxwell  equation.  Due  to 
the  electrochemical  reactions,  the  concentrations  of  both  re¬ 
actants  drop  to  zero  in  the  catalyst  layers.  Since  the  reactions 
occur  only  in  the  inlet  portion  of  the  catalyst  layers  in  this 
case,  the  thickness  of  both  catalyst  layers  may  be  reduced  to 
save  the  precious  catalyst  metal. 

The  distribution  of  water  vapor  concentration  in  terms 
of  relative  humidity  (RH)  is  presented  in  Fig.  3(b). 


Fig.  3.  Profiles  of  (a)  hydrogen  and  oxygen  concentrations,  c h-,  and  co2,  (b)  relative  humidity,  RH,  (c)  water  vapor  flux,  Afw>g,  (d)  liquid  water  flux,  7VWfi,  (e) 
temperature  difference,  AT  and  (f)  potential,  0,  across  the  cell  thickness  at  three  current  densities.  The  shaded  regions  denote  the  catalyst  layers,  and  the  results 
correspond  to  the  parameter  combination  given  in  Table  2. 
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Fig.  4.  Profiles  of  (a)  hydrogen  and  oxygen  concentrations,  ch2  and  co2,  (b)  relative  humidity,  RH,  (c)  water  vapor  flux,  AfW)g,  (d)  liquid  water  flux,  iVw>i,  (e) 
temperature  difference,  AT  and  (f)  potential,  0,  across  the  cell  thickness  at  three  current  densities.  The  shaded  regions  denote  the  catalyst  layers,  and  the  results 
correspond  to  pa  =  1  atm,  and  the  other  parameters  are  given  in  Tables  1  and  2. 


Recall  that  the  total  pressure  in  the  electrodes  is  assumed 
to  be  constant;  therefore,  the  concentration  of  the  water  va¬ 
por  increases  towards  the  catalyst  layer/electrode  interfaces 
to  balance  the  corresponding  decreases  in  reactant  concentra¬ 
tions  in  Fig.  3(a).  For  higher  current  density,  the  water  vapor 
concentration  near  the  catalyst  layer/electrode  interfaces  be¬ 
comes  higher,  which  is  instrumental  for  membrane  hydration. 
The  values  of  relative  humidity  are  greater  than  1  in  Fig.  3(b), 
indicating  that  water  condensation  occurs  in  both  electrode 
regions. 

Fig.  3(c)  and  (d)  shows  the  flux  distribution  of  water  va¬ 
por  and  liquid  water,  respectively,  across  the  cell.  Because 
of  the  above-mentioned  condensation  process,  /Vw.g  (Awj) 
decreases  (increases)  monotonically  with  x  in  the  electrode 
regions.  Note  that  water  vapor  is  present  only  in  the  electrode 
regions,  while  liquid  water  exists  throughout  the  cell.  In  the 
membrane  region,  an  increase  in  current  density  leads  to  in¬ 
crease  in  osmotic  drag,  which,  in  turn,  results  in  increase  in 
the  liquid  water  flux.  Rapid  increase  of  the  liquid  water  flux 


with  increasing  x  is  observed  in  the  cathode  catalyst  layer, 
owing  to  the  electrochemical  reaction  that  generates  liquid 
water  in  the  region.  Since  the  water  fluxes  are  assumed  to  be 
proportional  to  the  corresponding  reactant  mixture  flux  [13], 
both  Vw,g  and  Nw  \  increase  with  increasing  current  density 
in  the  electrode  regions.  The  same  assumption  is  responsi¬ 
ble  for  the  discontinuity  of  Nw j  at  the  interface  of  catalyst 
layer/anode  electrode,  and  a  more  accurate  model  account¬ 
ing  for  the  two-phase  flow  in  the  porous  media  needs  to  be 
incorporated  in  a  future  work. 

The  temperature  profiles,  T  —  Tc,  across  the  entire  cell 
are  shown  in  Fig.  3(e)  for  the  base  case  conditions  given 
in  Table  2.  Linear  variations  are  observed  in  the  anode  and 
cathode  electrodes,  indicating  that  heat  transport  in  the  re¬ 
gion  is  predominantly  through  conduction  with  negligible 
heat  generation  from  resistive  heating.  It  is  seen  that  a  small 
temperature  rise  (<0.2  K)  is  achieved  when  /  =  0.3  A  cm-2. 
With  increasing  current  density,  relatively  larger  tempera¬ 
ture  rise  (about  1.1  K  for  the  case  of  I  =  1.3  Acm-2)  and 
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highly  nonlinear  profiles  are  observed  in  the  membrane  re¬ 
gion,  owing  to  the  Joule  heating  from  the  ionic  current.  Since 
multiple  cells  are  stacked  together  in  typical  application,  the 
seemingly  small  temperature  rise  within  a  single  cell  and  the 
temperature  jump  across  the  cell  interface  caused  by  con¬ 
tact  resistance  may  accumulate,  thereby  causing  a  problem 
of  thermal  management. 

Fig.  3(f)  shows  the  electrical  potential  profiles  along  the 
cell  thickness.  A  value  of  1 . 19  V  is  set  as  the  reference  poten¬ 
tial  at  the  inlet  of  the  anode  electrode,  i.e.,  point  a  in  Fig.  1. 
The  overpotential  in  the  electrodes  is  negligible,  owing  to  the 
high  electrical  conductivity  of  the  electrode  material.  Signif¬ 
icant  activation  loss  and  ohmic  loss  are  seen  in  the  catalyst 
layers  and  the  membrane  region,  respectively.  It  must  be  men¬ 
tioned  that  the  anode  activation  loss  is  much  larger  compared 
to  the  results  in  the  literature  (e.g.,  [21  ]),  which  may  be  caused 
by  a  relatively  slow  hydrogen  electrochemical  reaction  rate 
adopted  in  this  study.  Note  that  uncertainties  exist  for  the 
value  of  the  reaction  rates  since  the  value  of  some  model  pa¬ 
rameters  (e.g.  Cjef  in  Eq.  (14))  are  not  well  documented  in 
the  literature.  Furthermore,  the  polarizations  increase  with 
increasing  current  density,  resulting  in  lower  cell  voltage  as 
shown  in  Fig.  2(a)-(d). 

Following  the  presentation  format  in  Fig.  3,  Fig.  4(a)- 
(f)  illustrates  the  effects  of  CO  concentration  on  the  dis¬ 
tribution  of  the  variables  of  interest  at  a  current  density 
I  —  0.5  A  cm-2,  and  the  parameter  combination  given  in 
Tables  2  and  1.  As  seen  previously  in  Fig.  3(a),  the  con¬ 
centration  gradient  of  FE  in  the  anode  electrode  in  Fig.  4(a) 
is  smaller  than  that  of  Ot,  due  to  the  larger  diffusivity  of 
H2  in  the  anode  gas  mixture.  A  small  decrease  of  ch2  in 
the  anode  electrode  is  observed  with  increasing  CO  concen¬ 
tration  since  the  total  pressure  remains  constant,  while  co2 
on  the  cathode  side  is  independent  of  the  CO  concentration 
on  the  anode  side.  In  Fig.  4(b),  the  relative  humidity  in  the 
anode  electrode  shows  small  decrease  with  increasing  CO 
concentration,  which  may  be  explained  by  the  reduction  of 
water  vapor  diffusivity  due  to  the  presence  of  the  CO  species. 
Again,  the  CO  concentration  in  the  anode  electrode  exerts  no 
influence  on  the  relative  humidity  in  the  cathode  electrode. 

The  water  vapor  flux  profile  in  Fig.  4(c)  remains  fixed  with 
respect  to  CO  concentration,  owing  to  the  previously  men¬ 
tioned  assumption  that  the  water  vapor  flux  is  proportional  to 
the  flux  of  the  reactant  mixture  (or  current  density).  Fig.  4(d) 
shows  that  liquid  water  flux  is  independent  of  the  CO  con¬ 
centration,  except  in  the  anode  catalyst  region.  A  decrease  in 
liquid  water  flux  in  the  anode  catalyst  region  is  observed  with 
increasing  CO  concentration,  which  is  due  to  the  increased 
consumption  of  liquid  water  during  the  electrochemical  reac¬ 
tion  denoted  by  Eq.  7.  Since  higher  CO  concentration  leads 
to  larger  concentration  and  activation  losses  in  the  catalyst 
region,  the  temperature  profile  in  Fig.  4(e)  shows  larger  tem¬ 
perature  difference  with  increasing  ppm  CO.  When  the  CO 
concentration  increases  from  10  to  25  ppm,  slight  increase 
of  T  —  Tc  is  seen,  while  relatively  large  increase  of  T  —  Tc 
occurs  as  ppm  CO  increases  further  to  100.  In  Fig.  4(f),  the 


overpotential  in  the  anode  catalyst  region  is  seen  to  increase 
with  the  increase  in  CO  concentration,  and  large  voltage  loss 
is  observed  when  ppm  CO  equal  to  100.  Note  that  the  ohmic 
loss  in  the  membrane  is  invariant  with  respect  to  the  CO  con¬ 
centration. 

With  the  overview  of  the  transport  and  electrochemical 
processes  illustrated  in  Figs.  3  and  4,  the  model  is  further 
used  in  a  parametric  study  to  evaluate  the  performance  of  the 
fuel  cell  over  a  wide  range  of  design  and  operation  condi¬ 
tions.  In  this  study,  four  quantities,  namely  the  cell  potential, 
Eceli-  the  power  density,  P,  the  maximum  temperature  rise, 
A7niax,  and  the  minimum  hydration  7,lnm  are  adopted  to  eval¬ 
uate  the  overall  fuel  cell  performance.  It  must  be  mentioned 
that  the  minimum  hydration,  is  determined  by  the  water 
vapor  activity  at  the  membrane/anode  catalyst  interface  [13]. 
Systematic  parametric  runs  are  performed  to  calculate  the 
four  performance  quantities  as  functions  of  the  parameters 
involved  in  the  model,  and  relevant  results  that  are  instru¬ 
mental  to  design  and  operating  optimization  are  presented  in 
the  remainder  of  this  section.  The  parameters  selected  in  the 
presentation  are  classified  into  two  categories:  (I)  the  operat¬ 
ing  parameters — current  density,  7,  CO  concentration  in  ppm, 
cell  temperature,  Tc,  anode  pressure,  pd,  and  relative  humid¬ 
ity  at  anode,  RHa  and  (II)  the  design  parameters — platinum 
loading,  ;«c,pt,  membrane  thickness,  fme,  and  porosity  of  the 
electrodes,  ee.  The  range  of  values  of  the  operating  and  design 
parameters  in  the  parametric  study  are  chosen  based  on  those 
commonly  used  in  the  literature  [1,8,9,13,21].  The  results 
from  the  parametric  study  are  used  for  developing  operating 
and  design  windows  based  on  the  maximization  of  the  power 
density  subject  to  practical  constraints. 

Fig.  5(a)— (d)  presents  the  performance  quantities  as  func¬ 
tions  of  the  current  density,  /,  for  different  values  of  CO  con¬ 
centration.  In  Fig.  5(a),  the  cell  potential,  £ceu,  decreases 
monotonic  ally  with  increasing  I  for  fixed  ppm  CO  and  with 
increasing  ppm  CO  for  fixed  current  density,  as  previously 
seen  in  Fig.  2.  When  the  ppm  CO  >  50,  the  limiting  current 
density  is  observed  to  be  reduced  to  around  0.6  A  cm-2,  and 
the  corresponding  cell  potential  reduces  to  a  low  value  of 
0.2  V.  Fig.  5(b)  shows  the  power  density,  P,  as  a  function  of 
the  current  density,  /,  and  ppm  CO.  Since  the  power  density 
P  is  defined  as  the  product  of  the  cell  potential  and  current 
density,  P  initially  increases  with  increasing  7,  for  fixed  CO 
concentration  in  Fig.  5(b).  With  further  increase  in  the  current 
density,  the  power  density  may  decrease  due  to  smaller  cell 
potential,  leading  to  the  formation  of  peaks  for  the  cases  of 
ppm  CO  <  50.  When  ppm  CO  >  100,  the  drop  of  the  cell  po¬ 
tential  is  not  significant  in  the  vicinity  of  the  limiting  current 
density,  and  no  power  density  peak  is  observed.  For  a  fixed 
value  of  current  density,  P  decreases  monotonically  with  in¬ 
creasing  CO  concentration,  owing  to  smaller  cell  potential 
caused  by  the  CO  poisoning  effect. 

The  maximum  temperature  difference  in  the  cell,  A7max, 
and  the  minimum  hydration  in  the  membrane,  Amm,  are  shown 
as  functions  of  7  and  CO  ppm  in  Fig.  5(c)  and  (d),  respec¬ 
tively.  The  maximum  temperature  difference  is  seen  to  in- 
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Fig.  5.  Parametric  effects  of  /  and  ppm  CO  on  (a)  cell  potential,  (b)  power  density,  (c)  maximum  temperature  difference  and  (d)  minimum  hydration.  The 
default  values  of  the  other  parameters  are  given  in  Tables  1  and  2. 


crease  monotonically  with  increasing  current  density,  owing 
to  more  Joule  heating  in  the  membrane  region  at  larger  I. 
For  fixed  value  of  current  density,  the  increase  in  CO  con¬ 
centration  leads  to  higher  overpotential  in  the  catalyst  region, 
which,  in  turn,  is  responsible  for  the  increase  in  Armax.  Since 
the  membrane  hydration  depends  primarily  on  the  concentra¬ 
tion  of  the  water  species,  Amin  is  seen  to  remain  a  constant 
value  of  14  for  the  range  of  current  density  and  CO  con¬ 
centration  considered  in  Fig.  5(d).  The  membrane  material 
is  in  a  well-hydrated  state  with  high  ionic  conductivity  as 
^■min  A  14. 

Following  the  presentation  format  in  Fig.  5,  Figs.  6-8  il¬ 
lustrate  the  effects  of  the  other  operating  parameters  on  the 
performance  quantities.  In  Fig.  6(a),  a  range  of  operating 
temperature,  Tc,  between  293  and  363  K  is  chosen  for  typical 
PEM  fuel  cell  operation.  For  fixed  value  of  Tc,  the  cell  po¬ 
tential  Eceii  decreases  monotonically  with  increasing  current 
density,  as  seen  previously  in  Figs.  2-5.  The  cell  potential  in¬ 
creases  when  the  operating  temperature  increases  from  293  to 
353  K,  which  may  be  explained  by  the  reduced  losses  in  the 
cell.  The  ionic  conductivity  increases  with  increasing  tem¬ 
perature,  leading  to  smaller  resistive  loss  in  the  membrane. 
The  transport  losses  are  also  reduced  since  the  diffusivities  in¬ 
creases  with  increasing  Tc .  Furthermore,  the  activation  loss  in 
the  catalyst  layer  is  smaller  at  higher  operating  temperatures 
[13].  However,  further  increase  in  operating  temperature  may 
cause  increased  partial  water  vapor  pressure,  leading  to  en¬ 


hanced  mass  transport  losses.  Consequently,  the  cell  voltage 
is  seen  to  slightly  decrease  as  Tc  increases  to  363  K. 

Fig.  6(b)  presents  the  power  density,  P,  as  a  function  of 
the  current  density,  /,  and  the  operating  temperature,  Tc.  The 
effect  of  the  operating  temperature  on  the  power  density  is 
similar  to  that  on  the  cell  potential,  i.e.,  P  increases  initially 
with  increasing  Tc,  and  starts  to  decrease  when  T  =  363  K. 
The  maximum  temperature  difference  in  Fig.  6(c)  decreases 
as  Tc  increases  from  293  to  353  K,  owing  to  the  decreased 
losses  mentioned  above;  a  slight  increase  in  A  7nlax  is  caused 
by  the  increased  transport  losses  when  Tc  =  363  K.  The  mini¬ 
mum  hydration  in  the  membrane,  Amin,  is  shown  as  a  function 
of  I  for  different  values  of  Tc  in  Fig.  6(d).  For  fixed  value  of 
Tc,  Amin  is  seen  to  decrease  monotonically  with  increasing  /, 
which  may  be  attributed  to  the  loss  of  the  water  species  at  the 
anode  catalyst/membrane  interface  caused  by  the  increased 
osmotic  drag  effect.  The  decrease  in  /,IIlm  is  more  pronounced 
at  lower  temperatures,  since  water  species  cannot  be  effec¬ 
tively  transferred  to  the  membrane  due  to  reduced  diffusiv¬ 
ities  with  temperature.  The  minimum  hydration  is  seen  to 
increase  monotonically  with  the  operating  temperature  for 
fixed  values  of  /,  owing  to  increased  water  vapor  pressure  at 
higher  temperatures. 

The  effects  of  the  anode  pressure  on  the  performance  quan¬ 
tities  are  shown  in  Fig.  7(a)-(d).  A  non-monotonic  trend 
is  again  observed  in  Fig.  7(a)  for  the  relation  between  the 
cell  potential,  £Ceii-  and  the  anode  pressure,  pa.  When  pa 
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Fig.  6.  Parametric  effects  of  /  and  Tc  on  (a)  cell  potential,  (b)  power  density,  (c)  maximum  temperature  difference  and  (d)  minimum  hydration.  The  default 
values  of  the  other  parameters  are  given  in  Tables  1  and  2. 


Fig.  7.  Parametric  effects  of  /  and  pa  on  (a)  cell  potential,  (b)  power  density,  (c)  maximum  temperature  difference  and  (d)  minimum  hydration.  The  default 
values  of  the  other  parameters  are  given  in  Tables  1  and  2. 
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(b)  Current  Density,  I  [A/cm2] 


Fig.  8.  Parametric  effects  of  I  and  RHa  on  (a)  cell  potential,  (b)  power  density,  (c)  maximum  temperature  difference  and  (d)  minimum  hydration.  The  default 
values  of  the  other  parameters  are  given  in  Tables  1  and  2. 


increases  from  1  to  3  atm,  higher  partial  pressure  (or  con¬ 
centration)  of  hydrogen  is  present  at  the  anode  catalyst  layer, 
leading  to  reduced  concentration  loss  or  higher  cell  potential. 
Further  increase  in  pa  results  in  reduced  water  vapor  diffu- 
sivity,  which  causes  lower  membrane  hydration.  The  higher 
ohmic  loss  in  the  membrane,  in  turn,  contributes  to  the  lower 
cell  potential.  Overall,  the  influence  of  pa  on  the  cell  poten¬ 
tial  is  not  as  significant  as  that  of  the  operating  temperature 
in  Fig.  6(a),  for  the  range  of  parameters  considered  in  this 
study. 

Fig.  7(b)  shows  a  similar  dependence  of  the  power  den¬ 
sity  on  the  anode  pressure  pd,  i.e.,  a  maximal  value  of  P 
is  achieved  at  a  moderate  pressure.  The  maximum  tempera¬ 
ture  difference,  A7nlax,  is  shown  as  a  function  of  the  current 
density,  7,  at  different  values  of  anode  pressure  in  Fig.  7(c). 
With  increasing  anode  pressure  p.d,  the  potential  losses  first 
decrease,  leading  to  the  decrease  in  A  Tmdx  since  less  heat  is 
generated  from  the  losses.  Further  increase  in  pd  results  in  in¬ 
creased  losses,  which  corresponds  to  the  increase  in  A  7nlax  at 
relatively  larger  pa.  The  variation  of  temperature,  however,  is 
observed  to  be  within  0. 1  K  for  the  range  of  pd .  The  minimum 
hydration,  kmin,  presented  in  Fig.  7(d)  shows  a  decreasing 
trend  with  increasing  anode  pressure,  owing  to  the  reduced 
water  vapor  diffusivity.  The  dehydration  of  the  membrane  is 
apparent  at  higher  anode  pressure  and  current  density,  and  a 
range  of  /,nnn  from  14  to  10  is  observed  in  Fig.  7(d). 


Since  PEM  fuel  cells  operated  under  conditions  of  low 
temperature,  high  current  density  and  high  pressure  are  prone 
to  membrane  dehydration  at  the  anode  side,  the  anode  feed 
stream  is  always  humidified  with  water  vapor.  Fig.  8(a)-(d) 
presents  the  effect  of  the  relative  humidity  in  the  anode  feed 
stream,  RHa,  on  the  performance  quantities  of  the  fuel  cell. 
In  Fig.  8(a),  the  cell  potential  increases  monotonically  with 
increasing  RHa,  on  account  of  improved  hydration  and  re¬ 
duced  ohmic  losses.  No  apparent  improvement  in  cell  poten¬ 
tial  is  observed  when  RHa  increases  from  1.0  to  1.1,  which 
may  be  explained  by  the  fact  that  the  supersaturated  water 
vapor  (RHa  =1.1)  will  quickly  be  condensed  in  the  anode 
electrode,  and  the  actual  water  vapor  concentration  in  the 
anode  catalyst  layer  is  equivalent  to  the  case  of  RHa  =  1 .0. 
Similarly,  the  power  density  increases  monotonically  with 
increasing  RHa,  as  shown  in  Fig.  8(b).  The  maximum  tem¬ 
perature  difference  in  Fig.  8(c)  decreases  with  increasing  an¬ 
ode  relative  humidity,  on  account  of  increased  hydration  and 
reduced  Joule  heating  in  the  membrane  region.  Relatively 
small  variation  of  the  temperature  (<0.2  K)  is  observed  for 
the  wide  range  of  RHa  considered.  The  minimum  hydration, 
7. min j  presented  in  Fig.  8(d)  shows  monotonic  increase  with 
increasing  RHa,  as  physically  expected.  However,  no  change 
of  Amjn  is  observed  with  increasing  7  for  a  fixed  value  of  RHa . 

It  is  evident  from  the  parametric  effects  presented  in 
Figs.  6-8  that  increasing  the  operating  temperature  can  gen- 
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erally  reduce  ohmic,  transport  and  activation  losses.  How¬ 
ever,  excessive  increase  in  Tc  leads  to  higher  water  vapor 
pressure,  which,  in  turn,  results  in  enhanced  transport  loss. 
Higher  anode  pressure  may  decrease  the  transport  loss,  how¬ 
ever,  membrane  dehydration  may  also  occur  due  to  decreased 
water  vapor  diffusivity.  Using  humidified  feed  stream  is  fa¬ 
vorable  to  the  membrane  hydration,  while  the  accompanying 
water  condensation  may  cause  flooding  of  the  electrodes.  An 
optimization  of  these  operating  parameters  is  therefore  nec¬ 
essary  to  determine  operating  regimes  that  carefully  balance 
these  competing  considerations.  In  this  study,  the  objective  of 
the  optimization  problem  is  considered  to  be  that  of  maximiz¬ 
ing  the  power  density,  subject  to  constraints  on  the  maximum 
temperature  difference,  the  minimum  membrane  hydration, 
and  the  maximum  cell  potential  as  explained  below.  The  goal 
of  the  optimization  problem  is  to  develop  the  ranges  of  fea¬ 
sible  current  density,  7,  as  function  of  the  other  operation 
parameters. 

The  optimization  scheme  considered  in  this  study  serves 
as  an  illustrative  example  to  introduce  the  methodology  for 
using  physics-based  models  for  cell  design  and  optimization. 
Based  on  other  operational  considerations,  alternative  opti¬ 
mization  problems  with  different  objective  function  and  con¬ 
straints  may  be  solved  following  the  methodology  outlined 
in  this  study.  The  roles  of  the  three  constraints  in  deciding 
the  bounds  on  the  operating  windows  are  as  follows: 

•  The  consideration  of  membrane  degradation  requires  that 
the  maximum  temperature  rise,  A  7max,  must  be  lower  than 
a  critical  value,  A.TCIit.  Since  A 7jnax  increases  with  in¬ 
creasing  current  density,  the  constraint  on  the  maximum 
temperature  corresponds  to  an  upper  bound  on  the  cur¬ 
rent  density.  It  must  be  pointed  out  that  the  present  one¬ 
dimensional  model  can  only  predict  the  across-the-cell 
temperature  rise,  which  is  significantly  lower  than  the  in¬ 
cell  temperature  rise  from  the  reactant  inlet  to  the  outlet. 
A  future  work  may  introduce  the  in-cell  temperature  rise 
constraint  by  adopting  a  two-  or  three-dimensional  fuel 
cell  model  in  the  optimization  study. 

•  The  membrane  hydration  is  calculated  as  the  ratio  of  the 
number  of  water  molecules  to  the  number  of  charge  sites 
in  the  membrane,  and  continuously  changes  from  0  for  a 
completely  dry  state  to  a  typical  value  of  16.8  for  a  fully 
hydrated  state  [10,11].  Since  the  ionic  conductivity  of  the 
membrane  material  increases  with  increasing  membrane 
hydration,  the  minimum  hydration  in  the  membrane,  /.mln, 
must  therefore  be  larger  than  a  lower  limit  /,cn[.  Since 
the  membrane  hydration  decreases  with  increasing  cur¬ 
rent  density,  the  constraint  on  the  minimum  membrane 
hydration  represents  a  second  upper  bound  on  the  current 
density.  Note  that  alternatively,  the  hydration  requirement 
may  be  posed  as  an  equality  constraint  on  /,mm  being  equal 
to  the  fully  hydrated  value.  The  results  presented  in  this 
study,  however,  are  based  on  the  inequality  constraint  on 


•  Since  the  current  density  decreases  with  increasing  cell 
potential,  the  size  and  the  capital  cost  of  the  fuel  cell  sys¬ 
tem  increase  for  the  higher  cell  potential.  The  cell  poten¬ 
tial  must  therefore  have  a  upper  bound,  Ttcrit,  determined 
based  on  the  capital  cost  consideration.  Low  cell  potential 
will  reduce  the  efficiency,  leading  to  higher  operating  cost. 
Thus,  the  cell  potential  must  also  have  a  lower  bound  deter¬ 
mined  from  the  efficiency  consideration.  The  cell  potential 
decreases  with  increasing  7,  consequently,  the  constraint 
on  the  upper  bound  of  cell  potential  forms  a  lower  bound 
on  the  current  density.  The  constraint  on  the  lower  bound 
of  cell  potential  merely  provides  an  additional  upper  bound 
on  the  current  density,  and  is  not  included  in  the  present 
optimization  scheme  for  the  purpose  of  conciseness. 

The  region  bounded  between  the  lower  bound  and  the 
smaller  of  the  two  upper  bounds  on  the  current  density  con¬ 
stitutes  the  operating  window  for  the  fuel  cell  based  on  the 
three  constraints.  The  optimum  current  density  variation  is 
one  that  lies  within  the  operating  window  and  maximizes  the 
power  density. 

The  optimization  problem  may  be  written  mathematically 


as 

Maximize  P  (21) 

subject  to: 

ATmdx  —  A  /'cut  <  0  (22) 

An  in  7-crit  A  0  (23) 

EceU  —  7?crit  <  0  (24) 


where  the  critical  values  in  this  study  are  chosen  as:  A7’crit  = 
1.0K,  Acr;t  =  14  and  Lcnt  =  0.75  V. 

Based  on  the  parametric  study  shown  in  Fig.  6,  the  con¬ 
struction  of  an  example  operating  window  on  the  current  den¬ 
sity,  7,  is  illustrated  in  Fig.  9(a)  as  a  function  of  the  operating 
temperature,  Tc .  Corresponding  to  the  constraint  on  the  max¬ 
imum  allowable  temperature,  shown  by  the  thick  straight  line 
in  Fig.  6(c),  the  upper  limits  of  the  current  density  for  differ¬ 
ent  values  of  Tc  are  indicated  by  the  solid  line  AB  in  Fig.  9(a). 
Similarly,  the  long-dashed  line  CD  and  the  chain-dashed  line 
EF  in  Fig.  9(a)  are  determined  based  on  the  constraints  on  the 
minimum  membrane  hydration  and  the  maximum  cell  po¬ 
tential,  shown  in  Fig.  6(d)  and  (a),  respectively.  The  shaded 
region  in  Fig.  9(a),  enclosed  among  the  three  constraints, 
is  identified  as  the  operating  window.  Values  of  the  current 
density  that  lie  inside  the  operating  window  will  ensure  that 
the  maximum  temperature  difference,  the  minimum  mem¬ 
brane  hydration  and  the  maximum  cell  potential  are  within 
the  prescribed  limits  that  correlate  to  desired  overall  cell  per¬ 
formance. 

While  satisfying  the  specified  constraints  is  a  necessary 
condition,  maximizing  the  power  density  simultaneously  is 
imperative  for  affordable  fuel  cells  with  small  size.  The  vari¬ 
ation  of  the  peak  value  of  the  power  density  with  Tc  (shown 
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(a)  Operating  Temperature,  Tc  [K] 


(c)  Anode  Relative  Humidity,  RH^ 


Fig.  9.  Operating  windows  as  functions  of  (a)  operating  temperature,  (b) 
anode  pressure  and  (c)  relative  humidity. 

by  the  thick  solid  curve  in  Fig.  6(b))  is  plotted  as  the  short- 
dashed  line  GH  in  Fig.  9(a).  This  line  represents  the  uncon¬ 
strained  current  density  as  a  function  of  operating  temper¬ 
ature,  Tc,  required  for  maximizing  the  power  density.  The 
superposition  of  the  unconstrained  solution  on  the  process¬ 
ing  window  is  used  to  identify  the  optimal  variation  of  the 
current  density  as  a  function  of  the  operating  temperature  so 
as  to  maximize  the  power  density,  P.  It  is  seen  in  Fig.  9(a) 
that  in  the  region  GI,  since  the  unconstrained  current  density 
values  for  maximizing  the  power  density  fall  above  the  upper 
bound  corresponding  to  the  Anml  constraint,  the  upper  bound 
constitutes  the  constrained  optimal  solution.  In  the  segment 
between  I  and  H,  the  unconstrained  optimum  current  density 
profile  passes  within  the  feasible  window,  and  as  a  result,  sat¬ 


isfies  all  the  specified  constraints.  The  constrained  optimum 
current  density  variation  with  the  operating  temperature  is  in¬ 
dicated  in  the  plot  by  the  thick  solid  line.  Since  the  constraints 
on  /.rnm  and  Zsmax  cannot  be  satisfied  simultaneously  when  Tc 
is  less  than  333  K,  neither  a  feasible  operating  window  nor  an 
optimal  solution  exist  for  Tc  <  333  K.  The  constrained  solu¬ 
tion  in  Fig.  9(a)  indicates  that  operating  conditions  for  fuel 
cell  systems  in  practice  mostly  correspond  to  off-peak  values 
of  power  density. 

Following  the  foregoing  procedure,  operating  windows 
and  optimal  solution  of  the  current  density  can  be  obtained 
for  the  other  operating  parameters  as  well,  and  the  results 
for  the  anode  pressure  and  the  anode  relative  humidity  are 
presented  in  Fig.  9(b)  and  (c),  respectively.  In  Fig.  9(b),  the 
optimal  solution  follows  the  unconstrained  optimum  current 
density  profile  for  pd  <  3.4 atm.  In  the  interval  3.4 atm  < 
Pa  <  5.0  atm,  the  unconstrained  optimum  resides  above  the 
upper  bound  corresponding  to  the  /,mm  constraint,  and  the  up¬ 
per  bound  constitutes  the  constrained  optimal  solution.  When 
pa  >  5  atm,  no  current  density  value  exists  that  satisfies  both 
the  ,A.min  and  Emax  constraints,  which  suggests  that  the  fuel 
cell  should  not  be  operated  with  anode  pressure  higher  than 
5.0  atm  for  the  given  parameter  combination.  Corresponding 
to  the  decrease  in  the  optimum  current  density,  the  power 
density  (not  shown  in  Fig.  9(b))  also  decreases  significantly 
with  increasing  anode  pressure,  which  may  be  explained  by 
the  fact  that  higher  pressure  reduces  the  species  diffusivities, 
leading  to  larger  concentration  losses.  The  decrease  in  power 
density  at  higher  anode  pressures  was  also  reported  by  Rowe 
and  Li  [13], 

In  Fig.  9(c),  the  optimal  solution  exists  when  the  relative 
humidity  is  greater  than  0.9.  When  RHa  <  0.9,  the  membrane 
hydration  may  become  small  to  violate  the  /,mm  constraint, 
or  the  cell  potential  may  be  high  enough  to  violate  the  £ceu 
constraint,  leading  to  poor  performance  of  the  cell.  The  results 
indicate  that  the  anode  feed  stream  needs  to  be  sufficiently 
hydrated,  as  physically  expected. 

The  results  presented  so  far  pertain  to  the  effects  of  the  op¬ 
erating  parameters  on  the  performance  of  the  fuel  cell,  which 
are  instrumental  to  the  derivation  of  optimum  operating  con¬ 
ditions.  It  is  further  illustrative  to  examine  the  effects  of  the 
design  parameters  to  gain  insight  into  the  optimum  design  of 
the  fuel  cell,  and  the  corresponding  results  are  presented  in 
Figs.  10-13.  Platinum  catalyst  loading  is  a  critical  parameter 
that  governs  the  available  sites  for  the  electrochemical  reac¬ 
tions.  The  dependence  of  the  polarization  curve  on  the  cath¬ 
ode  platinum  loading,  ;«c,pt,  is  illustrated  in  Fig.  10.  When 
mcpt  increases  from  0.32  to  3.00  mg  cm-2,  the  cell  poten¬ 
tial  is  observed  to  have  relatively  large  improvement  (about 
0. 1  V  for  all  the  current  density  values  considered),  owing 
to  the  increased  reaction  area.  Further  increase  in  mcpt  to 
8.00  mg  cm-2  only  yields  weak  increase  in  the  cell  potential, 
indicating  that  an  optimum  platinum  loading  may  be  deter¬ 
mined  by  simultaneously  evaluating  the  increases  in  perfor¬ 
mance  and  cost.  The  power  density  in  Fig.  10(b)  shows  a  sim¬ 
ilar  increasing  pattern  with  increasing  wc,pt,  as  physically  ex- 
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(b)  Current  Density,  I  [A/cm2] 
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Fig.  10.  Parametric  effects  of  1  and  /Hc,pt  on  (a)  cell  potential,  (b)  power  density,  (c)  maximum  temperature  difference  and  (d)  minimum  hydration.  The  default 
values  of  the  other  parameters  are  given  in  Tables  1  and  2. 
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Fig.  11.  Parametric  effects  of  I  and  fme  on  (a)  cell  potential,  (b)  power  density,  (c)  maximum  temperature  difference  and  (d)  minimum  hydration.  The  default 
values  of  the  other  parameters  are  given  in  Tables  1  and  2. 
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Fig.  12.  Parametric  effects  of  /  and  ee  on  (a)  cell  potential,  (b)  power  density, 
values  of  the  other  parameters  are  given  in  Tables  1  and  2. 


(c)  maximum  temperature  difference  and  (d)  minimum  hydration.  The  default 


pected.  With  increasing  platinum  loading,  the  activation  loss 
decreases,  leading  to  monotonic  decrease  in  the  maximum 
temperature  difference,  as  shown  in  Fig.  10(c).  Since  the  hy¬ 
dration  is  governed  by  the  water  species,  /nc,pt  is  observed  to 
have  no  obvious  effect  on 

Fig.  ll(a)-(d)  presents  the  effect  of  the  membrane 
thickness,  fme,  on  the  four  performance  quantities.  The 
increase  in  the  membrane  thickness  results  in  the  decrease 
in  cell  potential  in  Fig.  11(a),  owing  to  the  increased 
membrane  resistivity.  Similarly,  the  power  density  decreases 
monotonically  with  increasing  fme ,  as  shown  in  Fig.  11(b). 
With  the  increase  in  rme,  the  Joule  heating  in  the  membrane 
region  increases  monotonically,  resulting  in  the  increased 
A  7max  in  Fig.  11(c).  For  the  parameter  combinations  in 
Fig.  11(d),  the  minimum  hydration,  Amm,  is  observed  to  be 
independent  of  fme.  The  thicknesses  of  the  electrode  and 
catalyst  layers  show  similar  effect  as  that  of  fme,  and  the 
corresponding  results  are  omitted  here  for  brevity. 

The  effect  of  the  electrode  porosity,  et,  on  the  polarization 
curve  is  shown  in  Fig.  12(a).  Due  to  smaller  mass  transport 
loss  at  larger  electrode  porosity,  the  cell  potential  is  seen 
to  increase  monotonically  with  increasing  et.  For  the  cases 
with  high  current  density  (/  >  0.9  A  cm-2),  the  mass  trans¬ 
port  loss  constitutes  a  large  portion  of  the  total  loss;  thus, 
the  increase  in  ee  leads  to  relatively  large  increase  in  cell  po¬ 
tential,  as  seen  in  Fig.  12(a).  Similarly,  the  power  density  in 


Fig.  12(b)  shows  monotonic  increase  with  increasing  <?e,  and 
again,  larger  increases  are  observed  at  large  values  of  cur¬ 
rent  density.  The  maximum  temperature  difference,  Armax, 
increases  as  the  electrode  porosity  increases,  which  may  be 
attributed  to  the  increasing  Joule  heating  in  the  high-porosity 
electrode  region  (Fig.  12(c)).  The  minimum  hydration  in 
Fig.  12(d)  remains  fairly  independent  of  the  electrode  poros¬ 
ity  and  the  current  density.  However,  for  the  case  ee  =  0.6, 
a  slight  decrease  in  /.rnm  is  observed  with  increasing  current 
density,  I. 

Again,  following  the  procedure  in  Fig.  9,  operating 
windows  and  optimal  solution  of  the  current  density  are 
obtained  for  the  design  parameters,  mCiPt,  hue  and  ee  as 
depicted  in  Fig.  13(a)-(c),  respectively.  In  Fig.  13(a),  the 
optimum  solution  line  follows  the  unconstrained  optimum 
current  density  profile  for  mc,pt  <  1.6  mg  cm-2.  When 
mc,pt  >  1.6  mg  cm-2,  the  unconstrained  optimum  solution 
exceeds  the  upper  bound  corresponding  to  the  A  7max 
constraint,  and  the  upper  bound  denotes  the  constrained 
optimal  solution.  For  the  entire  range  of  platinum  loading 
considered,  the  optimum  solution  of  current  density  exists 
and  only  slightly  varies  between  1.15  and  1.25  Acm-2.  The 
unconstrained  optimal  solution  line  in  Fig.  13(b)  is  observed 
to  be  close  to  the  two  upper  bounds,  and  the  optimum 
current  density  has  relatively  small  variation  (between 
1.15  and  1.30  A  cm-2)  with  respect  to  the  membrane 
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<c)  Electrode  Porosity,  e 
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Fig.  13.  Operating  windows  as  functions  of  (a)  cathode  platinum  loading, 
(b)  membrane  thickness  and  (c)  electrode  porosity. 


thickness.  However,  the  optimum  power  density  (not  shown 
in  Fig.  13(b))  exhibits  significant  increase  with  decreasing 
membrane  thickness  (see  Fig.  11(b)),  which  indicates  that 
thinner  membrane  is  preferred  in  practical  applications.  The 
operating  window  for  the  electrode  porosity  is  presented  in 
Fig.  13(c),  where  the  optimum  solution  line  follows  the  un¬ 
constrained  optimum  current  density  profile  as  et  <  0.415. 
The  optimum  solution  is  given  by  the  maximum  temperature 
difference  constraint  in  the  interval  0.415  <  ee  <  0.562,  and 
by  the  minimum  hydration  constraint  when  et  >  0.562. 

The  results  in  this  section  provide  a  comprehensive  anal¬ 
ysis  of  a  PEM  fuel  cell  operation,  and  illustrate  a  methodol¬ 
ogy  for  deriving  operating  envelopes  and  optimum  process 
designs  using  physics-based  simulation  of  the  cell.  Although 


the  operating  windows  presented  are  specific  to  the  constraint 
values  chosen,  other  constraints  may  be  incorporated  follow¬ 
ing  the  approach  outlined.  Note  that  the  optimum  solutions  in 
Figs.  9  and  13  are  obtained  by  changing  one  of  the  operating 
or  design  parameters  while  fixing  the  values  of  the  remain¬ 
ing  ones.  These  specific  solutions  determine  the  values  of  the 
current  density  that  maximize  the  power  density  for  the  par¬ 
ticular  parameter  combination  chosen  for  the  study.  Optimum 
parameter  combinations  will  be  reported  in  a  future  work  by 
coupling  the  present  model  with  numerical  optimization  pro¬ 
cedure  [23]. 


4.  Conclusions 

A  comprehensive  one-dimensional  steady-state  model  for 
a  single  fuel  cell  operated  with  reformate  fuel  is  presented, 
and  a  systematic  parametric  investigation  is  conducted  over 
a  wide  range  of  the  fuel  cell  operating  and  design  conditions. 
Using  the  process  simulation  results,  optimum  current  den¬ 
sity  is  obtained  for  the  first  time  as  a  function  of  various  op¬ 
erating  and  design  parameters.  For  the  specific  optimization 
scheme  and  parameter  combinations  adopted,  it  is  found  that 
higher  operating  temperature  Tc  and  anode  relative  humidity 
RHa,  and  lower  anode  pressure  pa  are  required  to  maximize 
the  power  density.  The  optimum  current  density  does  not 
have  significant  change  with  the  two  design  parameters  con¬ 
sidered  in  the  study,  i.e.,  the  platinum  loading,  /«cp[,  and  the 
membrane  thickness,  fme.  With  the  increase  in  the  electrode 
porosity,  the  temperature  difference  and  the  membrane  hydra¬ 
tion  may  become  limiting  factors  to  the  cell  power  density. 
The  results  provide  direct  usable  information  for  optimum 
operating  and  design  of  PEM  fuel  cells.  The  methodology 
for  determining  the  operating  windows  and  optimum  design 
configurations  may  be  readily  applied  to  other  constraints  and 
cell  design  specifications. 
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